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We analyze a mean-field model of electrons on a square lattice with two types 
of interaction: forward scattering favoring a d-wave Pomeranchuk instability and a 
BCS pairing interaction driving d-w&ve superconductivity. Tuning the interaction 
parameters a rich variety of phase diagrams is obtained. If the BCS interaction 
is not too strong, Fermi surface symmetry breaking is stabilized around van Hove 
filling, and coexists with superconductivity at low temperatures. For pure forward 
scattering Fermi surface symmetry breaking occurs typically via a first order transi- 
tion at low temperatures. The presence of superconductivity reduces the first order 
character of this transition and, if strong enough, can turn it into a continuous one. 
This gives rise to a quantum critical point within the superconducting phase. The 
superconducting gap tends to suppress Fermi surface symmetry breaking. For a rel- 
atively strong BCS interaction, Fermi surface symmetry breaking can be limited to 
intermediate temperatures, or can be suppressed completely by pairing. 
PACS: 71.10.Fd, 71.10.Pm, 74.20.Mn 

I. INTRODUCTION 

Usually the Fermi surface of an interacting electron systems respects the point-group 
symmetry of the underlying crystal lattice. However, electron-electron interactions may 
also lead to Fermi surface deformations which break the orientational symmetry sponta- 
neously. From a Fermi liquid viewpoint this can happen via a "Pomeranchuk instability" , 
that is, when Pomeranchuk's stability condition^ for the forward scattering interactions 
is violated. 

Interactions favoring a symmetry-breaking Fermi surface deformation with a ci-wave 
order parameter, where the surface expands along the fc^-axis and shrinks along the k y - 
axis (or vice versa), are present in the t-J,— Hubbard,- 1 ^ and extended Hubbard 5 model on 
a square lattice. These models therefore exhibit enhanced "nematic" correlations, which 
also appear in the context of fluctuating stripe order. 6 Signatures for such correlations 
have been observed in various cuprate superconductors.- In particular, they provide a 
natural explanation for the relatively strong in-plane anisotropy observed in the magnetic 
excitation spectrum of YBa 2 Cu30 y .-^ 

Fermi surface symmetry breaking competes with superconductivity. In the t-J model 
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the rf-wave Fermi surface deformation instability is overwhelmed by (i-wave pairing. This 
is indicated by slave-boson mean-field theory^ and has been confirmed recently by a vari- 
ational Monte Carlo calculation. 10 However, enhanced nematic correlations remain.— The 
competition of superconductivity and Fermi surface symmetry breaking is more delicate 
in the two-dimensional Hubbard model. Renormalization group calculations in the sym- 
metric phase suggest that the superconducting instability is always stronger than the 
Pomeranchuk instability— but these calculations do not exclude the possibility of coex- 
istence of the two competing order parameters in the symmetry-broken phase. Indeed, 
coexistence of <i-wave superconductivity and <i-wave Fermi surface symmetry breaking has 
been obtained near van Hove filling from a weak coupling perturbation expansion for the 
symmetry-broken ground state of the Hubbard model.— 

To elucidate the interplay and competition of Fermi surface symmetry breaking and 
superconductivity in a more general setting, and to classify possible scenarios, we analyze 
in the present work a mean-field model allowing for both instabilities with a tunable 
strength for each. The model describes itinerant electrons on a square lattice with two 
types of interaction: a reduced BCS interaction driving <i-wave superconductivity and a 
purely forward scattering interaction driving <i-wave Fermi surface symmetry breaking. 

The properties of the mean-field model without BCS interaction, where the electrons 
interact only via forward scattering ("f- model"), have been clarified already earlier.— 
The main results can be summarized as follows. Fermi surface symmetry-breaking occurs 
below a transition temperature T c which forms a dome-shaped line as a function of the 
chemical potential /i, with a maximal T c near van Hove filling.— 1 ^ The phase transition 
is usually first order at the edges of the transition line, and always second order around 
its center.— £LML The <i-wave compressibility of the Fermi surface is however strongly en- 
hanced even near the first order transition down to zero temperature.— Adding a uniform 
repulsion to the forward scattering interaction, the two tricritical points at the ends of the 
second order transition line are shifted to lower temperatures. For a favorable choice of 
hopping and interaction parameters one of the first order edges can be replaced completely 
by a second order transition line, leading to a quantum critical point—' Fluctuations at 
and near the quantum critical point destroy fermionic quasi-particle excitations, leading 
to non-Fermi liquid behavior.—^ 

Adding an attractive <i-wave BCS interaction to the f-model leads to a variety of 
qualitatively distinct phase diagrams, depending on the interaction strength. If the BCS 
interaction is not too strong, Fermi surface symmetry breaking is stabilized around van 
Hove filling, and coexists with superconductivity at low temperatures. In the presence 
of a pairing gap it is easier to realize Fermi surface symmetry breaking via a continuous 
phase transition at low temperatures than without. In particular, a quantum critical point 
connecting superconducting phases with and without Fermi surface symmetry breaking 
at zero temperature is obtained for a suitable choice of interactions. For a relatively 
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strong BCS interaction, Fermi surface symmetry breaking can be limited to intermediate 
temperatures, or can be suppressed completely by pairing. 

The article is structured as follows. In Sec. II we introduce the mean-field model and 
describe the self-consistency equations for the order parameters. The phase diagrams and 
other results are presented in Sec. III. A conclusion follows in Sec. IV. 

II. MEAN-FIELD MODEL 

We analyze itinerant electrons on a square lattice interacting via forward scattering 
and a reduced BCS interaction, described by a Hamiltonian of the form 

H = J2&* + H l + H $i (!) 

k 

where nk = counts the spin-summed number of electrons with momentum k. The 

kinetic energy is due to hopping between nearest and next-nearest neighbors on a square 
lattice, leading to the bare dispersion relation 

e k = — 2t(cos k x + cos k y ) — At' cos k x cos k y . (2) 

The forward scattering interaction reads 

H i = 2Z ^ kk ' ™ krik ' ' ( 3 ) 

k,k' 

where L is the number of lattice sites, and the function /kk' has the form 

/kk' =u- g/dkdv , (4) 

with coupling constants u > and gj > 0, and a function with d x 2_ y 2-wave symmetry 
such as <ik = cosk x — cosk y . This ansatz mimics the structure of the effective interaction 
in the forward scattering channel as obtained for the t-J— and Hubbard 3 model. The 
uniform term originates directly from the repulsion between electrons and suppresses the 
(uniform) electronic compressibility of the system. The c?-wave term enhances the d- 
wave compressibility and drives spontaneous Fermi surface symmetry breaking. In the 
Hubbard model it is generated by fluctuations, while in the t-J model the nearest neighbor 
interaction contributes directly to a <i-wave attraction in the forward scattering channel. 
The BCS interaction has the form 

Hj = j Vkk ' c kT c -ki c -k'|Ck'T • (5) 

k,k' 

It is a reduced interaction in the sense that it contributes only in the Cooper channel, 
that is, when the total momentum of the interacting particles vanishes. For the matrix 
element Vkk' we choose a separable <i-wave attraction 

Vkk' = ~9c dkd k > (6) 
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with g c > 0, which corresponds to the dominant term in the Cooper channel for the 
two-dimensional Hubbard and t-J model. 

Inserting n k = (n k ) + Sn^ into if/, and c^ci^ = (c k ^.c_ k |) + S(c^c_ k ^) into Hf, and 
neglecting terms quadratic in the fluctuations, one obtains the mean-field Hamiltonian 

Hmf = [ eknk + ( AkC k T c -k| + h - c -) ~ "y W ~ Ak ( c k T c -k|) • ( 7 ) 

k 

Here e k = e k + £e k is a renormalized dispersion relation, which is shifted with respect to 
the bare dispersion by 5e k = L^ 1 ^ k , /kk'(^k') = un + r/d^, where n = L^ 1 X^k( n k) is 
the average particle density, and 

V = -j;J2 d ^) (8) 

k 

is our order parameter for Fermi surface symmetry breaking. It vanishes as long as the 
momentum distribution function (n k ) respects the symmetry of the square lattice. The 
superconducting gap function is given by A k = j Vkk' (c_ k /jC k /-|-) = Ad k , where 

A = -|E^ c -^t>- (9) 

k 

For the reduced interactions if/ and Hj the mean-field decoupling is exact in the ther- 
modynamic limit. Feynman diagrams describing contributions beyond mean- field theory 
have zero measure for L — > oo. 

The mean-field Hamiltonian is quadratic in the Fermi operators and can be diagonal- 
ized by a Bogoliubov transformation. For the grand canonical potential per lattice site, 
uj = L^Q, we obtain 



(77, A) = - A £ l g[2 cosh(/?£ k /2)] + i + !^! + 1i „-^-^, (10) 



where (3 is the inverse temperature, = (£ k + 1 A^ 2 ) 1 / 2 , and £ k = e k — /i. The stationarity 
conditions du/di] = and du/dA = yield the self-consistency equations for the order 
parameters 



k k 

and 



- ^y4^tanh^, (12) 
L ^ 2£ k 2 ' v ; 

k k 

respectively. The condition du/dn = (at fixed /i) yields the equation determining the 
density 

n = l-iy| k -tanh^. (13) 
k K 
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III. RESULTS 

We now show results obtained from a numerical solution of the mean-field equations. 
For the ratio of hopping amlitudes we choose t'/t = —1/6. The bare dispersion e£ has 
saddle points at k = (n,0), (0,7r), leading to a logarithmic van Hove singularity in the 
bare density of states at e = —4t' = —2t/3. All the results presented in the figures are 
for u = (no uniform contribution to forward scattering), but we will discuss the effects 
of a finite u in the text. In the following we set t — 1, that is, all results with dimension 
of energy are in units of t. 

In Fig. 1 we show the transition temperature Tf(fi) for Fermi surface symmetry break- 
ing in the absence of pairing (A = 0) for gf — 1, and the critical temperature for super- 
conductivity T c (n) in the absence of Fermi surface symmetry breaking (rj = 0) for various 



choices of g c . As discussed in detail in Refs. |l5| and[l6|, a symmetry-broken Fermi surface 



is stabilized below a dome-shaped transition line, with a maximal transition temperature 
near van Hove filling. The transition is first order at the edges of the transition line and 
second order around its center. The critical temperature for superconductivity is also 
maximal near van Hove filling, but T c (/i) remains finite for any fi (as long as the band is 
partially filled), and the transition is always of second order. Near van Hove filling the 
transition temperatures Tf and T c are of the same order of magnitude for g c slightly above 
gf — 1. Note that in the weak coupling limit one would obtain Tf <C T c for comparable 
gj and g c , since log(l/7}) oc gj 1 for gf — > 0, 16 while log(l/T c ) oc gc 1 ^ 2 for g c — > at van 
Hove filling, due to the logarithmic divergence of the density of states, and the additional 
logarithm in the Cooper channel, as can be seen from the gap equation (12). 

We now discuss results for the full mean-field model, allowing also for coexistence of 
the two order parameters. In Fig. 2 we show the low temperature region of the phase 
diagram in the (fi, T)-plane for gj = 1 and a relatively weak BCS coupling, g c = 0.7. 
Fermi surface symmetry breaking suppresses T c and remains essentially unaffected by 
the (relativeley small) superconducting gap. The suppression of T c occurs since Fermi 
surface symmetry breaking splits the van Hove singularity, reducing thus the density of 
states at the Fermi level. However, superconductivity cannot be eliminated completely, 
since a logarithmic Cooper singularity survives for any reflection invariant Fermi surface. 
The phase diagram thus exhibits three types of first order transitions between phases 
with a symmetric and a symmetry-broken Fermi surface: between two normal states, 
between a superconductor and a normal state, and between two superconducting states. 
Continuing the (then metastable) phase with a symmetric Fermi surface beyond the first 
order transition line leads to a diverging <i-wave compressibility at the ficticious second 
order transition line "T| nd " also shown in the plot. 

For larger g c the energy scale for superconductivity (gap and T c ) increases, and effects of 
the superconducting gap on Fermi surface symmetry breaking become more pronounced, 
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see Fig. 3 (here g c = 0.9). In particular the first order transition line Tf(fi) is shifted 
toward the center of the symmetry-broken region, and approaches the ficticious second 
order line "T| nd " . In Fig. 4 we show the A-dependence of the "reduced" Landau energy 
uj(A) = w[r/ min (A), A] — u[r] min (0), 0] for two points in the phase diagram which are close 
to each other, but on opposite sides of the first order transition between a superconducting 
state with a symmetric Fermi surface and a normal state with Fermi surface symmetry 
breaking; 77 min (A) minimizes lu(t],A) for fixed A. Note that ?7 m i n is zero for large A; the 
kink in lo(A) is due to the discontinuous onset of r] for small A. The /^-dependence of the 
order parameters rj and A is shown for various temperatures in Fig. 5. The jump of rj 
at the first order transition induces a counter jump of A. For high T, superconductivity 
is suppressed completely by Fermi surface symmetry breaking (Fig. 5c), while for lower 
temperatures coexistence of the order parameters A and r] is realized (Figs. 5a and 5b). 
The temperature dependence of the order parameters is shown for /i = —0.7 (near van 
Hove filling) in Fig. 6. The increasing superconducting gap A leads to a decrease of rj upon 
lowering the temperature below T c . Superconductivity smears the single particle states 
over an energy range of order A, and thus suppresses the energy gain from a Fermi surface 
deformation. Fermi surface symmetry breaking is thus suppressed by the superconducting 
gap. 

Although the system is not critical at the first order transition from a symmetric to 
a symmetry-broken Fermi surface, it is close to criticality in the sense that the (i-wave 
compressibility Kd is strongly enhanced by the forward scattering interaction. For the case 



of pure forward scattering (f-model) this was shown already in Ref. [16|. In the presence 



of superconductivity with a gap function Ak, the ci-wave compressibility is given by 

«d = i ~ > ( 14 ) 

1 - 9f Krf 



where 

H,l — — / u-k 



(15) 



2E* ( cosh f«k) 2 El 

is the <i-wave compressibility in the superconducting state in the absence of forward scat- 
tering (gf = 0). The enhancement of due to g$ is thus given by the "Stoner factor" 
S = (1 — gj k®)^ 1 . In Fig. 7 we plot the inverse Stoner factor S^ 1 along the right first 
order transition line (approached from the symmetric phase) up to the tricritical temper- 
ature Tj T1 for various choices of g c . It becomes clear that S is enhanced significantly by 
superconductivity at low temperatures. In particular, for g c = 0.9 the system is very close 
to criticality. 

Enhancing g c beyond g c = 0.9, the first order transition lines between the states with 
symmetric and symmetry-broken Fermi surfaces are successively replaced by a continuous 
phase transition. In Fig. 8 we show the phase diagram for g c — 1. Here Fermi surface 
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symmetry breaking occurs via a continuous transition at the lowest temperatures, well 
below T c . In particular, there is a continuous quantum phase transition at T = 0. The 
first order lines are connected to continuous transition lines both at the high and low 
temperature ends. The low temperature ends are tricritical points, where the quadratic 
and the quartic coefficient of the reduced Landau energy u)(rj) = U)[r), A min (^)] both vanish. 
By contrast, at the high temperature ends the quartic coefficient of uj(rf) jumps from a 
negative to a positive value. This discontinuity is due to the onset of A below T c . Note 
that the high temperature ends are close to the tricritical points found for smaller g c , 
such that a small jump of the quartic coefficient can turn its sign. For g c = 1.12 the first 
order transition has disappeared completely from the phase diagram (see Fig. 9), and the 
transition between symmetric and symmetry-broken Fermi surfaces is always continuous. 
The transition lines for Fermi surface symmetry breaking and superconductivity intersect 
in tetracritical points, where both quadratic coefficients of u(j],A) vanish. Enhancing 
g c further leads to a progressive suppression of Fermi surface symmetry breaking, in 
particular at lower temperatures, where the superconducting gap is getting large. For 
g c = 1.2, Fermi surface symmetry breaking is eliminated completely by superconductivity 
at low T, while it still survives in a small region at intermediate temperatures, see Fig. 10. 
For even larger g c the region with a symmetry-broken Fermi surface shrinks further until 
it disappears completely from the phase diagram. 

Adding a uniform contribution u > to /kk'j Eq. (4), leads to a suppression of first 
order transitions into a phase with a symmetry-broken Fermi surface, making thus con- 
tinuous transitions easier. This trend was already observed and explained in detail for the 
case of pure forward scattering.— For small g c , the tricritical points are shifted to lower 
temperatures by a finite u, and the first order transition line moves closer to the ficticious 
second order transition. The gradual replacement of the first order line by a second order 
for increasing g c is accelerated for u > 0. For example, for gf = 1, u — 10, and g c = 0.9 
the phase diagram looks qualitatively as the one in Fig. 9, with Fermi surface symmetry 
breaking always occuring via a continuous transition. 

The (effective) interaction resulting from the Hubbard or t-J model contains also an 
s-wave component in the Cooper channel. In case of coexistence of superconductivity 
with a (i-wave Fermi surface deformation, this leads to a small s-wave contribution to the 
gap function Ak, in addition to the dominant d-wave term.-^ 

IV. CONCLUSIONS 

We have solved a mean-field model for itinerant electrons moving on a square lattice 
with two types of interactions: an interaction in the forward scattering channel favoring 
a <i-wave shaped symmetry-breaking Fermi surface deformation and a reduced BCS inter- 
action with (i-wave symmetry. Making different choices for the interaction parameters, a 
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rich variety of possible phase diagrams has been found. 

For pure forward scattering Fermi surface symmetry breaking occurs typically via a first 
order transition at low temperatures.— 1 ^ The presence of superconductivity reduces the 
first order character of this transition and, if strong enough, can turn it into a continuous 
one. This gives rise to a quantum critical point within the superconducting phase. The 
superconducting gap tends to suppress Fermi surface symmetry breaking. For a certain 
choice of parameters one finds reentrant behavior, where Fermi surface symmetry breaking 
is stabilized at intermediate temperatures, while it is suppressed by the pairing gap at 
low temperatures. If superconductivity is too strong, Fermi surface symmetry breaking 
disappears completely from the phase diagram. 

In microscopic models the relative strength of forward scattering and pairing interac- 
tions is determined by the microscopic interactions. In the t-J model slave-boson mean- 
field^ and variational Monte Carlo^ calculations show that pairing prevents Fermi surface 
symmetry breaking, but there are strongly enhanced correlations indicating that the model 
is close to a (i-wave Pomeranchuk instability. 11 This corresponds to the case of a relatively 
large g c in our phenomenological mean-field model. In the weakly interacting Hubbard 
model coexistence of superconductivity and Fermi surface symmetry breaking has been 
found around van Hove filling at T = within second order perturbation theory— The 
available numerical results indicate that Fermi surface symmetry breaking occurs via a 
continuous transition in this case, as in the phase diagrams in Figs. 8 or 9. 

It would clearly be interesting to analyze how order parameter fluctuations modify the 
mean-field results. A renormalization group calculation by Vojta et ai— suggests that a 
quantum critical point for orientational symmetry breaking in a ci-wave superconductor 
is destabilized by fluctuations, leading possibly to a first order transition. 

Acknowledgments 

We thank Marijana Kircan for a critical reading of the manuscript and for valuable 
comments. 



1 I. J. Pomeranchuk, Sov. Phys. JETP 8, 361 (1958). 

2 H. Yamase and H. Kohno, J. Phys. Soc. Jpn. 69, 332 (2000); 69, 2151 (2000). 

3 C. J. Halboth and W. Metzner, Phys. Rev. Lett. 85, 5162 (2000). 

4 I. Grote, E. Kording, and F. Wegner, J. Low Temp. Phys. 126, 1385 (2002); V. Hankevych, 
I. Grote, and F. Wegner, Phys. Rev. B 66, 094516 (2002). 

5 B. Valenzuela and M. A. H. Vozmediano, Phys. Rev. B 63, 153103 (2001). 

6 S. A. Kivelson, E. Fradkin, and V. J. Emery, Nature 393, 550 (1998). 



9 



7 S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, J. M. Tranquada, A. Kapitulnik, 
and C. Howald, Rev. Mod. Phys. 75, 1201 (2003). 

8 V. Hinkov, S. Pailhes, P. Bourges, Y. Sidis, A. Ivanov, A. Kulakov, C. T. Lin, D. P. Chen, 
C. Bernhard, and B. Keimer, Nature 430, 650 (2004). 

9 H. Yamase and W. Metzner, Phys. Rev. B 73, 214517 (2006). 

B. Edegger, V. N. Muthukumar, and C. Gros, Phys. Rev. B 74, 165109 (2006). 

1 H. Yamase, Phys. Rev. Lett. 93, 266404 (2004). 

2 C. Honerkamp, M. Salmhofer, and T. M. Rice, Eur. Phys. J. B 27, 127 (2002). 

3 A. Neumayr and W. Metzner, Phys. Rev. B 67, 035112 (2003). 

4 H.-Y. Kee, E. H. Kim, and C.-H. Chung, Phys. Rev. B 68, 245109 (2003). 

5 I. Khavkine, C.-H. Chung, V. Oganesyan, and H.-Y. Kee, Phys. Rev. B 70, 155110 (2004). 

6 H. Yamase, V. Oganesyan, and W. Metzner, Phys. Rev. B 72, 035114 (2005). 

7 W. Metzner, D. Rohe, and S. Andergassen, Phys. Rev. Lett. 91, 066402 (2003). 

8 L. Dell' Anna and W. Metzner, Phys. Rev. B 73, 045127 (2006); |corrd-mat/061172~3 

9 M. Vojta, Y. Zhang, and S. Sachdev, Phys. Rev. Lett. 85, 4940 (2000); Int. J. Mod. Phys. B 
14, 3719 (2000). 



10 



0.25: 
0.2 
0.15; 

0.1 ; 

0.05 : 




-2 



-1.5 



-1 -0.5 





T-2nd 

i f I 




f . 


• 


T-itri 

h 


■ 


T c (g c =\A) 




K(g =1.2) 




r c (&=i.o) 


A 


Tc(g c =0.7) 



0.5 



FIG. 1: Critical temperature Tf(fi) for Fermi surface symmetry breaking in the absence of 
superconductivity (A = 0) for gf = 1 and critical temperature for superconductivity T c (/x) 
in the absence of Fermi surface symmetry breaking (77 = 0) for various choices of g c . The 
superconducting transition is always of second order. Fermi surface symmetry breaking occurs 
via a first order transition at temperatures below the tricritical points Tj T1 , and via a second 
order transition above. 
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FIG. 2: Phase diagram in the (/_t, T)-plane for gf = 1 and g c = 0.7. The symbols 77 and A indicate 
which order parameters are finite in the various regions confined by the transition temperatures. 
Fermi surface symmetry breaking occurs via a first order transition in the temperature range 
shown in the plot. Continuing the (then metastable) phase with a symmetric Fermi surface 
beyond the first order transition line leads to a diverging d-wave compressibility at the ficticious 
second order transition line "T? nd " also shown in the plot. 




11 



0.12 




-0.8 -0.6 



FIG. 3: Phase diagram in the ( / u,T)-plane for gf = 1 and g c = 0.9. The first order transition 
line for Fermi surface symmetry breaking is very close to the ficticious second order line in this 
case. 
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FIG. 4: Reduced Landau energy a>(A) = u;[ry m i n (A), A] — uo[r] m i n (0), 0], minimized with respect 
to r], as a function of |A|. The interaction parameters are 5/ = 1 and g c = 0.9 as in Fig. 3, the 
temperature is T = 0.09. The two choices of \i correspond to two points close to but on opposite 
sides of the first order transition line between a superconducting state with a symmetric Fermi 
surface and a normal state with a d-wave deformed Fermi surface. 
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FIG. 5: Order parameters rj and A as a function of \i for various temperatures, T = 0.05, 0.07, 
0.09. The interaction parameters are gf = 1 and <? c = 0.9 as in Fig. 3. 




FIG. 6: Temperature dependence of the order parameters rj and A for fi = —0.7 and interaction 
parameters as in Fig. 3. 
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FIG. 7: Temperature dependence of the inverse Stoner enhancement of the <i-wave compressibil- 
ity along the right first order transition line between the phases with symmetric and symmetry- 
broken Fermi surface. Interaction parameters are gf = 1 and g c = 0, 0.7, 0.9. The tricritical 
point at the high temperature end of the transition line (T| n = 0.130) is the same in all cases. 
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FIG. 8: Phase diagram in the (fi, T)-plane for gf = 1 and 5f c = 1. The ficticious second order 
transition line "T? nd " is so close to the first order line Tl st that it is hidden by the latter. 
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FIG. 10: Phase diagram in the (/x, T)-plane for gf = 1 and g c = 1.2. All transitions are of second 
order. 



